The magnetic susceptibility of the t-J model at low hole doping 
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We compute the dynamical magnetic susceptibility of the t-J model in its commensurate flux phase 
at low hole doping. We compare the calculations with experiments and exact diagonalization studies. 



It has long been recognized that an essential ingredient 
in any theory of high T c superconductivity is a descrip- 
tion of the accompanying magnetic behavior jjj. There- 
fore, a great deal of effort has been invested in study- 
ing the magnetic behavior of the t-J model, one of the 
simplest models thought to capture the physics of the 
high temperature superconductors S. Because the t-J 
model is very strongly correlated, it is not easily solvable. 
Nonetheless, two of the present authors applied a novel 
calculational technique to this model and the resulting 
optical response function agreed very well with exact di- 
agonalization studies ||,[|]- In this paper, we extend this 
calculational formalism in order to calculate zero tem- 
perature magnetic properties of the t-J model. We find 
quantitative agreement with exact diagonalization stud- 
ies at low dopings and show that our results are consis- 
tent with several experimental properties of the high T c 
materials. 

We will study the t-J Hamiltonian 
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where < i,j > denotes the Hermitian sum over nearest 
neighbor pairs, and no lattice site may be doubly occu- 
pied. In the gauge theory of the t-J model, the spin and 
charge degrees of freedom are treated as separate while 
the Gutzwiller constraint of having no doubly occupied 
sites is enforced by a gauge field degree of freedom. The 
corresponding Lagrangian is 
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a charge one spinless excitation at site i, and j\ a is a 
neutral operator that creates a spin a excitation at site 
i. The time component of the U(l) gauge field at site j 
is given by the Lagrange multiplier <pj while the phase of 
the Hubbard-Stratonivich variable \jk is the integral of 
the spatial component of this gauge field along the link 
joining sites j and k. 

Our calculational procedure requires us to define a sad- 
dle point of this gauge theory |H). We first define the 
mean-field Hamiltonians for the Fermi and Bose sectors 
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where we have assumed that Xjk 
independent magnitude \ Q given by 



1 



-— (0IW/ + 74I0) 



(5) 



where c\ a — fj a bi, b\ is a bosonic operator that creates 



It is important to note that the value of Xo is twice that 
predicted by the saddle point procedure. This bandwidth 
enhancement is a known effect of Gutzwiller projection 
The commensurate flux saddle point is defined by 
choosing mean- field values for Ojfc, the value of the gauge 
field along the link joining sites j and k, such that each 
plaquette of the lattice encloses (1 — S)/2 flux quanta, 
where S is the percentage of empty sites. We are able to 
investigate small fluctuations about these mean-field val- 
ues for the flux per plaquette by employing perturbation 
theory. One of the reasons we choose the commensurate 
flux saddle point is that it minimizes the total energy 
of these mean field Hamiltonians in the relevant doping 
(2) range P| . For the fermions, Hf is a Hofstadter Hamilto- 
nian [pTwith M/N = (1 — 5)/2 flux quanta per plaquette 
(M, N: relatively prime integers). In order to perform 
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calculations in the bose sector, a statistical transmuta- 
tion is employed that allows us to treat the bosons in 
their fermionic representation This assumption sta- 
bilizes the theory at zero temperature because the bosons 
are no longer degenerate and have a corresponding mean- 
field Hamiltonian with an energy gap. In addition, this 
approximation is consistent with known variational re- 
sults @. 

Having defined the underlying formalism, let us now 
calculate the spin correlation function 5(q, uf). We first 
calculate the mean field spinon polarization bubble, 



4H = ^E 
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is short-ranged and therefore immune to "screening" ef- 
fects. Therefore, U is not an adjustable parameter in this 
calculation. The vertex correction is easily performed, 
and the resulting expression for the corrected suscepti- 
bility is 



XqM = (1 - S) —f r . 

1 + E/XqM 
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where the doping-dependent prefactor is implicit in the 
relation 

S S = \fU<rUf ]P {b 3 b]) ~ (1 - S)\f] a [<rUfifi ■ (8) 



where 5+ 



Eq, and the sum is carried out over all excited states 
ofH f . 
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FIG. 1. Comparison of 5(q, w) calculated from Eq. (8) 
(solid lines) with the exact diagonalization results of ref. |12|] 
(dashed lines) at two different momenta for two, three and 
four holes on a 16-site cluster for J ft = 0.4. Both calculations 
use n — J/2. 

Next we perform a vertex correction in order to ac- 
count for the gauge field fluctuations that impose the 
Gutzwiller constraint. We approximate the gauge field 
interaction as an instantaneous onsite repulsive poten- 
tial of strength U. We require that U = 1.78 Jxo, the 
value found in variational studies of the Heisenberg an- 
tiferromagnet p| . This value of U, without Gutzwiller 
projection, generates the exact amount of antiferromag- 
netic order which minimizes the energy of the projected 
flux ground state. We assume that the repulsive core of 
the gauge field interaction between spinons is relatively 
unaffected at low doping. The reason is that the core 
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FIG. 2. Top: the spin structure factor S(q) for 5 — 3/47 
(solid line) compared with the one-hole exact diagonalization 
results of ref. jjj (o, J/t = 0.4) and ref. jjj (□, J/t = 0.25; 
+, J/t — 0.5). Inset: Brillouin zone, with special points la- 
beled. Bottom: 5(q) for 5 = 5/41 (solid line) compared with 
the two-hole results of ref. Q (o, J/t = 0.4) and ref. jljj 
(+, J/t — 0.4). Inset: S(q) calculated for (top to bottom) 
8 w 1/16, 1/8, 3/16, and 1/4. Note that our 5(q) is indepen- 
dent of t/J. 

Let us now compare our calculation of S'(q, uS) = 
-ImXq(w) with the exact diagonalization results of ref. 
|Q. Fig. [l] shows 5(q, w) at q = Q (left columns) and 
q = |Q (right) for dopings of 2/16, 3/16 and 4/16 (top, 
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middle and bottom respectively). The exact diagonal- 
ization values are shown as dotted lines. The solid lines 
were obtained from Eq. (|7]). We wish to point out the 
following basic similarities: 

• There is a single coherent peak with a characteristic 
energy scale of J at q = Q. The energy of this pole 
increases substantially with doping. 

• Elsewhere in the zone, (e.g. q = |Q) the spectral 
weight is not especially doping-dependent, nor is it 
very coherently distributed. 




5 q(n,Tl) 

FIG. 3. Left: The discommensuration qs in the lowest 
energy spin excitation (solid line) compared with the exper- 
imental data of ref. (TsJ . Due to an antiferromagnetic insta- 
bility below 8 — 0.02 (shaded), we have not calculated qs in 
this regime. Top right: contour plot of S(q, too) calculated 
for 8 = 15/107 ~ 0.14 (J/i=0.4, u = 1.14 J, r] = 0.02 J). 
Darker shading indicates a larger value of i5(q, Uo)- Bottom 
right: the Lai.seSro.i-iCuCU neutron scattering data of ref. 
JL 6[ . We have added the solid line as a guide to the eye. 
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FIG. 4. Solid line: The spin gap frequency u)q/J as 
a function of doping, for J/t — 0.4. Diamonds (o): The 
Neel temperature Tjv of the antiferromagnetic phase for 
La2-i Sr{Cu04, taken from the experimental data of ref. [ pc| , 
for dopings 8 < 0.02. The vanishing spin gap frequency as 
5 — ► 5j indicates the instability to the antiferromagnetic state 
(shaded) characterized by Tat. 



The unphysical sharp structure in both calculations 
has been smoothed out considerably through the use of a 
large value of 77 (= J/2). With the value of U held fixed, 
there are no adjustable parameters in our equation for 
S(q,oj). Furthermore, we have not scaled the y axes. 

We shall now calculate 5(q) = J Q duiS(q,uj) to see 
how the integrated spectral weight is distributed over the 
Brillouin zone. This gives an idea of the extent of mag- 
netic correlation at momentum q. Fig. ^compares S(q) 
with the exact diagonalization results of refs. ]l3|-|i~5|]. As 
expected, the correlation is greatest near Q and is sup- 
pressed elsewhere in the zone. While the exact diagonal- 
ization results are weakly dependent on J/t at Q, our 
calculation is independent of J/t. The doping depen- 
dence and absolute magnitude of S(q) are both repro- 
duced quite well. Again, we emphasize that there are no 
adjustable parameters. 

Let us now address the issue of discommensuration in 
S(q, lj). The plateau in 5(q) near Q widens upon doping, 
an effect clearly visible in the inset in the lower half of 
Fig. 0. The widening plateau is an indication of discom- 
mensuration of the antiferromagnetic correlations away 
from the ordering vector Q. The extent of this is usually 
characterized by a scalar quantity qs, where 

<W = -|<U-Q|, (9) 

and qd is a momentum at which the greatest magnetic 
scattering occurs. Neutron scattering experiments per- 
formed on L^-aSr^CuC^ have given values of qs at 
S = 0.075 and 0.14 (c.f. Fig. |, bottom right panel) 
Jl6[ . These values are shown in the left panel of Fig. |j, 
along with our calculation of qs (S) . The upper right panel 
shows how we estimate qs from S(q, to). This particular 
example shows a doping of 15/107 ~ 0.14. The lowest 
energy poles in 5(q, oj) occur at an energy uj = 1.14 J, 
and are distributed around Q in a ring-like manner (dark 
shading). At higher dopings, the ring becomes increas- 
ingly diamond-shaped, with more spectral weight along 
the faces than at the corners. To be consistent with the 
experiments, we take q^ in the (it, 0) direction and cal- 
culate qs as shown. 

While the magnitude of the discommensuration does 
agree with the experimentally obtained values from 
La2-aSr,5Cu04, the distribution does not. The experi- 
ments clearly show that discommensurate momenta ap- 
pear at four points in the (0, ±tt) and (±tt, 0) directions, 
not in a ring around Q. This disagreement between our 
calculation and the experimental data does not neces- 
sarily call into question the validity of our calculation. 
Exact diagonalization studies of the t-J Hamiltonian are 
as yet unable to determine which type of discommensu- 
ration it exhibits. Resolution of this issue will have to 
wait until clusters of large enough size can be studied. 
For now, we note that the discommensuration effect that 
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we calculated is due to the existence of a finite concen- 
tration of empty lattice sites. Because this distribution 
is uncorrelated (i.e. there is no charge order), all direc- 
tions look equivalent to spin excitations of momentum 
q. Therefore, a ring of discommcnsurate momenta forms 
around the ordering vector (n, tt). We expect this to 
change if the ground state is modified to include axis- 
aligned charge ordering, a modification which would re- 
quire the addition of another Hartree-Fock parameter. 
The formation of such "stripes" is thought by some to be 
an essential ingredient of high T c superconductivity [ p"7[ , 
and it has indeed been observed in experiments |l8[ |. If 
such ordering were found to occur along the crystalline 
axes it would cause discommensuration primarily along 
the (±7r, 0) and (0, ±7r) directions. 

It is well known that the cuprate superconductors be- 
come antiferromagnets in the limit of extreme underdop- 
ing. The doping at which this occurs has been estimated 
in several calculations |Ts|| . Our ground state becomes 
unstable to antifcrromagnctism at dopings below 2%, 
which happens to be exactly the dopant concentration 
below which long range antiferromagnetic order appears 
in La2_5Sr,5Cu04 p0| . The instability is indicated by 
a vanishing gap to spin excitations as <5 — > 8^ , where 
5 C ~ 0.02. Above 5 C , the lowest lying spin excitation re- 
quires a finite energy ujq. Fig. |] illustrates the boundary 
between the stable regime (6 > S c , ujo > 0) and the un- 
stable regime (S < S c , uj = 0). We have included the 
values of the Neel temperature measured in ref. [^0| for 
La 2 _aSr5Cu04 at dopings S < 0.02. 

While our calculation of the minimum doping of the 
spin disordered phase S c agrees with experimental stud- 
ies of La2_5Sr5Cu04, our calculation of the gap to spin 
excitations loq in this phase does not. Neutron scatter- 
ing investigations of the high T c materials indicate a gap 
to spin excitations with an energy scale closer to the su- 
perconducting gap than to J (glj. Both our calculation 
and the exact diagonalization work thus far greatly over- 
estimate the value of wo- While this may be a finite 
size effect in the exact diagonalizations, we attribute it 
to a overly large mean-field energy gap that is expected 
to be reduced by broadening and damping effects such as 
those considered in ref. |2|] with respect to the fractional 
statistics gas. 

To summarize, we have calculated the zero tempera- 
ture magnetic properties of the t-J Hamiltonian in its 
commensurate flux phase. We find quantitative agree- 
ment with exact diagonalization calculations of the spin 
correlation function ]l2] | and the spin structure factor 
-]l5|. We also find discommensuration of magnetic 
order away from (it, n) comparable in magnitude to that 
observed in La2_5Si'5Cu04 Jl6| , but the details remain 
to be clarified. The instability to an antiferromagnetic 
phase occurs at a value of 5 close to that found in ex- 
periments [^0| and other calculations [ ]l9fl . However, the 
gap to spin excitations in the commensurate flux phase 



is approximately five times larger than indicated by ex- 
periment |2~j] ]. Nevertheless, the t-J model seems to be 
well described by a gauge theory of new particles identi- 
fied with the spin and charge coordinates of the original 
electrons. We therefore present these calculations as a 
further test of the gauge theory of spinons and holons 
and as a justification for continued efforts to remedy its 
shortcomings. 
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